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Abstract 

The microscopic correlation functions of non-chiral random matrix models with complex eigen- 
values are analyzed for a wide class of non-Gaussian measures. In the large-iV limit of weak 
non-Hermiticity, where N is the size of the complex matrices, we can prove that all /c-point cor- 
relation functions including an arbitrary number of Dirac mass terms are universal close to the 
origin. To this aim we establish the universality of the asymptotics of orthogonal polynomials in 
the complex plane. The universality of the correlation functions then follows from that of the kernel 
of orthogonal polynomials and a mapping of massive to massless correlators. 
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1 Introduction 



The Universality of correlation functions is one of the crucial properties of random matrix models. 
In any of their applications, where we refer to |l]] for a review, the mapping between the underlying 
field theory and the matrix model as an effective model is only dictated by the global symmetries of 
the field theory. On the matrix model side the invariance under this symmetry, given for example by 
unitary transformations U(N), does not entirely fix the matrix model action and allows for a large 
variety of terms. The statement of universality is that in general the quantities of interest such as 
matrix eigenvalue correlation functions are the same for a certain class of matrix model actions in the 
large- N limit. The simplifying choice of for example a Gaussian therefore does not matter. 

One has to distinguish between different kinds of universality. In the macroscopic large- iV limit the 
eigenvalues of the random matrices which for example correspond to energy levels are kept unsealed. 
In this limit the oscillatory behavior of the correlation functions is smoothed and all two- and higher 
fc-point functions are universal, to all orders in an expansion in 1/N 2 Pj. Such universality results 
find applications in Quantum Gravity [§] , condensed matter physics of conductors @ or more recently 
in the correspondence between supersymmetric Yang-Mills theory and low-energy String Theory [H]. 

In this work we concentrate on a different universality called microscopic. Here, eigenvalues and 
correlation functions are rescaled with the size of the matrices N, magnifying a particular part of 
the energy spectrum. Microscopic universality states that correlation functions are the same for a 
class of matrix model actions when measured with respect to the same mean level spacing between 
the eigenvalues. The eigenvalue spectrum of random matrix models is typically compact, similar for 
example to lattice gauge theories where large momenta are cut off due to the finite lattice spacing. In 
the simplest case the matrix model spectral density given by the Wigner semi-circle on an interval. 
We distinguish between the origin, bulk and edge region of the spectrum where random matrix model 
universality was proven in ||, || and |7j respectively (for a review see ||). The origin of the spectrum 
is particularly interesting for applications to chiral symmetry breaking in low-energy Quantum Chro- 
modynamics (QCD) || because of the relation between the order parameter, the chiral condensate, 
and the spectral density of QCD Dirac operator eigenvalues at the origin 113]. Here, universality has 



been shown for a large class of polynomial [p and non-polynomial potentials pM, at criticality [12 



and in the presence of massive Dirac terms [113 1 . For classes with orthogonal and symplectic symmetry 
proofs including massless and massive Dirac terms have been given in [14] and [15| respectively. 

All we have said so far is valid for random matrices with real eigenvalues only. However, more 
recently applications for models with complex eigenvalues have occured as well. To give some examples 
we mention the fractional Quantum Hall Effect [16], two-dimensional charged plasmas [17, |18| |, QCD 
with chemical potential [19], as well as the inverse potential problem in integrable systems pp| . We 
note that the chiral ensembles also contain complex matrices || [| but remain in the class of real 
eigenvalue models. 

The number of different symmetry classes among complex matrix models is particularly rich [21]. 
Furthermore, one has to distinguish between the large- N limit of strongly and weakly non-Hermitian 
matrices. While in the former case results date back to Ginibre ]E2T| the latter limit has been defined 



rather recently [^] and correlation functions have been calculated in [24, 25, 2£, £7]]. The class 
of microscopic correlators in the weak non-Hermiticity limit plays an important role as it permits to 
interpolate between the correlations of real eigenvalues and those of complex eigenvalues at strong non- 
Hermiticity. Our studies are devoted to the microscopic origin scaling limit at weak non-Hermiticity, 
where up to date almost nothing was known about universality. The only exception is the microscopic 
spectral density which has been shown to be universal in [28] using supersymmetric techniques for 
ensembles of sparse matrices with independently distributed entries and finite second moment. While 
supersymmetric techniques become quite involved for higher correlators we will exploit the existence 
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of orthogonal polynomials in the complex plane [16, 27] and analyze their asymptotic properties. Due 



to their universality and the power of the method of orthogonal polynomials [29] we can derive the 
universality of all /c-point functions. Along the way of our analysis we will encounter similarities 
between the complex matrix model and the two-matrix model with two independent matrices as being 



studied for example in [ 30 1 . 



We will restrict ourselves to the simplest model, the complex extension of the Hermitian one-matrix 



model as introduced and solved in |16|, f24j. In 25] the same model including an arbitrary number 



of Dirac masses was solved in the limit of weak and strong[| non-Hermiticity and its relevance for 
three dimensional QCD with chemical potential was pointed out. In order to show also microscopic 



universality of the chiral model with complex eigenvalues introduced and solved very recently in [27], 
we would need to know higher order 1/N corrections to our method. The chiral model is relevant 
for QCD with chemical potential in four dimensions. Apart from its universality it would be very 
interesting to see if it is indeed in the same class of models as 19]. 



The article is organized as follows. In Section |2j we define our non-Gaussian complex matrix model 
together with all relevant quantities at finite-iV. In the following Section || a second order differen- 
tial equation for the orthogonal polynomials in the complex plane is derived in order to determine 
their asymptotic behavior. This leads to the universality of all fc-point correlation functions of com- 
plex eigenvalues in the microscopic origin scaling limit at weak non-Hermiticity in Section ||. Our 
conclusions and future prospects are presented in Section [5]. 

2 The model 

The partition function of the matrix model we wish to study is defined as 

z (2N f ) ^ / dJdjj Y[det[(J-im f )(J^ + im f )\ exp -NixV(J, J f ) , (2.1) 
f=i 

y(j,jt) ^ t^{ jj] - T 2 {j2 + j]2) ) + 2:Eff( j2fc + jt2fc )- ( 2 - 2 ) 



k=2 



The complex N x N matrix J is parameterized as 



J = H + \ l ^ — -iA, rG [0,1], (2.3) 

V 1 + T 

where H and A are Hermitian. The potential V(J, J') consists of a quadratic, Gaussian part depending 
on the parameter r as well as of an arbitrary but fixed number of higher even powers in J and of 
maximal degree 2d. The Gaussian part is fixed by requiring a Gaussian weight for both matrices H 
and A at equal variance (1 + r)/(2N). When enlarging the measure to study universality we have to 
make a choice which terms we allow. Instead of allowing for arbitrary even potentials for the matrices 
H and A respectively, inserting J and jt according to eq. ( |2.3| ), we have only introduced even powers 
of in J and J\ The reason why we have to exclude higher order powers of mixed terms, (JJ^) k - 2 , 
goes as follows. In order to write eq. fl2.ip in terms of eigenvalues we have to do the following unitary 
transformation, J = U(Z + R)U^ Here U is unitary, Z = diag(zi, . . . , zn) is the diagonal matrix 

1 In the massless limit the one- and two-point function were already calculated in |l7j at maximally strong non- 
Hermiticity, usins the picture of a charged plasma. The arbitrary charge at the origin there corresponds to massless 
Dirac terms in 1231. 
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of complex eigenvalues of J, and R is a strictly upper triangular matrix]^]. Since products of strictly 
upper triangular and diagonal matrices remain traceless, the independence drops out in all powers 
trJ 2fc and similarly in its Hermitian conjugate. In tr( Jjt) = tr(Z Z* -\- RBJ) the -R-dependence persists 
but can be integrated out because it is Gaussian. Adding higher powers tr(Jjt) fc mixed terms between 
Z, Z^ and RR^ would occur and the matrix R could no longer be integrated out. This is why we 
restrict ourselves to the form of the potential eq. (p.2|). A similar class of potentials has also been 
considered in |2Cj. In terms of complex eigenvalues %=i,. „,jv the partition function eq. ( ^l| ) thus reads 



r(2iV/) 



N 

n 



d 2 z 



V(z,z*) 



1 ( « 



Nf 

• n 

r 



im f \ 2 e - NV ^ 



i> \A N (Z!,...,Z N ) 



-(Z + Z 



* 2- 



k=2 



(2.4) 



(2.5) 



with Ajv(zi, • • • , zjsi) = Wkyiizk — z{) being the Vandermonde determinant. We have dropped the 
constants coming from the integration over the matrices U and R. The /c-point correlation functions 
among the eigenvalues Zj are defined in the usual way 2£j: 



R 



(2N f ) 



N 



(zi, 



■Zk) 



(N-k)l z ( 2N f) 



(2.6) 



N 



N 



d 2 z k +i 



Mil 

i=i \/=i 



1 2 -NV(zi,z*) 
imf\ e v 3 i' 



|Ajv(^i, 



>z N ) 



As it has already been pointed out in [15, 25] the correlation functions with and without the Nf quark 
masses are related. Absorbing the mass terms into the Vandermonde the following relation^] has been 
derived (2^]|: 



R 



(2N f ) 



N 



(Zl, ... ,Zfc) 



R 



{ N +Nf {zi, . . .,z k ,imi, . . .,im Nf ] 



(2.7) 



R N+Nf {imi, . . .,im Nf ) 

Consequently we only have to calculate the correlators at Nf = and prove their universality. The 
universality of the massive correlation functions then follows from eq. (2.7) at large-A. 

The fc-point correlation functions eq. (|2.6| ), massless or massive, can be obtained in the standard 



(2i 



way 2£] using orthogonal polynomials. These are defined in the complex plane as 

J d 2 z e- NV ^*) Pk ( z ) Pl ( z *) = S M , 



were from now on we suppress the superscript (2Nf = 0). The kernel of orthogonal polynomials is 
then given by 

N-l 



K N {z 1 ,z* 2 ) = e" 



1=0 

and the correlation functions follow from it |2T 



R N (zi,. . . ,z k ) 



Jet jK N (z u z*)] . 



(2.9) 



(2.10) 



2 We could also diagonalize J by a biunitary transformation, J = UAV^ . This would diagonalize products (JJ') k but 
not J 2 and 2 in the Gaussian part. 

3 At finite- JV the relation holds exactly when the iV-dependence in the weight is absorbed into the eigenvalues, e~ v( - z ' z ' . 
In eq. (2.21) of ref. [Ea] an additional N is missing in the subscript of the denominator. 
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Already at finite- iV the orthogonal polynomials eq. (2.8) possess a matrix or eigenvalue representation 
for an arbitrary potential V(J,J^) eq. (2.2). For that purpose we introduce a parameter dependent 
partition function with matrices of smaller size n x n, where n < N: 



Z n {t) = / (dJdjt) nxn exp ~\xV{J,rf) 



As one can easily convince oneself the following representation for orthogonal polynomials 
holds in the complex plane: 



(2.11) 
ill also 



Pn{z) 



1 



Zn(t 



n • 
N i 



(dJdJ^) nxn det(z-J) exp -NtrV(J,J r ) = ( det(z - J) ) Zn{t= n_) ■ (2.12) 



Here, they are normalized as = z n + . . . . In order to obtain orthonormal polynomials as in the 

definition (|2.8|) we divide by their norm, h n , to obtain 



P n {z) = h n ~ 2 P n (z) = h n ~ 2 (det(z- J)) 



(2.13) 



For the issue of universality we will also have to keep track of the (t = -^-dependence of the poly- 
nomials in the asymptotic large-n limit keeping t = jj finite. We repeat that eqs. ( |2.8| ) - ( p. 13 ) also 
hold in the presence of 2Nf masses in the measure. 

The resolvent for the partition function eq. ( |2.11| ) is defined as 



G n (z;t) = 



Z n {t) 



(dJdf) n 



1 

xn tr 

n z 



J 



exp 



n 



trF (J, J f ) 



-tr— 

n z 



J 



z n (t) 



It can also be expressed through the 1-point function, the spectral density, by 



G n (z;t) 



- I ,hr 

n 



R n (w) 



w 



(2.14) 



(2.15) 



This representation is an alternative to calculate in particular the macroscopic large-n spectral density 
from a saddle point analysis. We can deduce the following relation between density and resolvent in 
the complex plane which is different for real eigenvalues (see eq. ( |3.1l| ) below): 



1 



7T 



d z *G n (z;t) 



n 



-RJz) 



(2.16) 



3 The asymptotic of orthogonal polynomials 

In this section we will derive a differential equation for the orthogonal polynomials from the represen- 
tation eq. ( |2.13j ). The universality of the asymptotic polynomials and in consequence of the correlation 
functions then follow from it. The fact that we will take the large-n limit at weak non-Hermiticity 
will allow us to borrow some results from the proof in the real case Q despite the fact that the mi- 
croscopic correlations spread out into the complex plane in our case. The polynomials defined in eq. 
( p.S| ) have parity symmetry because of the even potential chosen. We will therefore aim at a second 
order differential equation that closes on polynomials of equal parity. 
We start by taking the derivative of P n {z) from eq.( 2.13j ) 



d z P n (z) 



-tr- 

n z 



1 



nG n ( z; t 



n 
N 



det(z - J) 
P n (z) + nhn 1 ( -lr- 



(3-1) 



n z 



J 



det{z - J) 
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where we have introduced the connected part (com) of the expectation value. At large-re expectation 
values are known to factorize and the second term will be suppressed by -=jr with respect to the first. 
We could in principle also determine the asymptotic of P n (z) from the above equation as it has been 
previously done in || in the real case. However, in order to obtain the correct behavior, sine and 
cosine as it will turn out, we would have to know also subleading corrections in N to fix their parity 
as in the real case |(|. Furthermore it is known that for orthogonal polynomials of real eigenvalues an 
exact, finite-re differential equation exists §§ which only closes on the polynomials as an equation of 
second order. Taking one more derivative we arrive at 



<):r n (z) - nil-G,, [z-^)P n {z) + n 2 G n [z-^) 2 P n {z) (3.2) 



^ v conn 1 / 1 1 \ coml 

+ nG n (z; hn 1 ( -tr -det(z-J)) + ( -tr -det(z-J)) 



./V/ \ Tt Z J I ( n \ \ Tt Z J I o- i n \ 

which is exact at finite n. At large-n the connected expectation values in the second line are sup- 
pressed by —j as follows from general counting arguments in matrix models with unitary invariance. 
Nevertheless the last term is of 0(1) since the derivative may produce another factor of n. We thus 
end up with the following second order differential equation: 



9 2 z Pn(z) 



2 n\ 2 ..,../ /; 



n z G n (z;-) + nd z G n (*; - ) + 0(1) 



Pn{z) . (3.3) 



We now proceed exactly as in the universality proofs for real eigenvalues. The coefficient functions in 
front of the polynomials in eq. ( |3.3| ) are smoothed and thus replaced by their values in the macroscopic 
large-n limit. The precise meaning of smoothing is defined as follows. The orthogonal polynomials 
P n (z) obey a so-called recursion relation where zP n -\{z) is reexpressed in terms of linear combinations 
of polynomials Pi(z) for I < n. On the real line a three step recursion relation holds for arbitrary 
polynomial potentials, which is no longer true in the complex plane. Here, the degree of the recursion 
explicitly depends on the degree 2d of the potential eq. (|2.5|) , similarly to the two- matrix model |3(| . 
Smoothing is defined by the assumption that the recursion coefficients r\ approach a single function 
in the variable t = 4?. Let us give an example. For the Gaussian potential given by V(z,z*) = 
1 _ 1 7 _ 2 (zz* — tj(z 2 + z* 2 )) the orthogonal polynomials are known to be Hermite polynomials in the 

complex plane [JO]], P n (z) = J ^j(l — t 2 )~^.t^ He n {z^J N/t), and they obey the following recursion 
relation^] 

zP^r(z) = ^P^ auss (z) + rP^r (*) • (3-4) 

Obviously the recursion coefficients r n = n/N approach a smooth limit, r n n±i,±2 . — ► r(t) = t. 
Taking this limit we obtain the following expression for the Gaussian resolvent outside the support a 
of the complex eigenvalues 

G G ^(z;t) = \[^~ , * e C\a , (3.5) 

which coincides with at t = 1. Here, we have used the fact that the resolvent eq. ( |2.15| ) can be 
expressed through the kernel via eq. ( |2.10| ) . Expanding the pole and smoothing the non- vanishing even 
powers, JcPze-^z^PtWPtiz*) = {Tt) k ( 2 *), where we can follow the calculation for real eigenvalues 



given e.g. in M, we arrive at eq. (3.5). It is the macroscopic resolvent outside the support as can be 



5t 



From eq. (3.4) one can see that the Christoffel-Darboux formula for the kernel eq. (f^Sl) no longer holds for r < 1 
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seen from the vanishing of eq. ( |2.16|) , d z *G Gauss (z;t) = 0. For our potential eq. (|2,5|) the recursion 
relation fl3.4| ) extends down to P n -2d+i(z)- We assume that the recursion coefficients have a smooth 
limit and that the procedure leads to the macroscopic resolvent, G(z, t), and its derivative d z G(z, t) in 
eq. ( |3.3[ ). Because of the complicated form of the recursion relation for our general potential we use 
another derivation of the macroscopic resolvent using loop equations [|j . Exploiting the invariance of 
the partition function eq. ( 2.11| ) under the change of variables J — > J + e/(z — J) and its Hermitian 



conjugate (H.c.) leads to the following equation, after requiring d e Z(t)\ ( 




n 



-tr 



(djV(J,J))- 



J 



0: 



+ H.c. 



(3.6) 



It contains mixed terms as J* ~ l /(z — J) which cannot be written solely in terms of eigenvalues, as 
explained in Section ^. This mixing can be partially removed for the Gaussian part of the potential 
with the help of the identity 







n 



-tr 



(djV(J,J^))- 



1 



Jt 



+ H.c. . 



(3.7) 



It follows from the invariance of Z(t) under J -> J + e/(z* - jt). Adding r times eq. ( |3.6| ) to eq. 
( O ) leads to 




1 d 

j+It 



2 ^ 

k=2 



2fc-l 



+ 



J 



+ H.c. 



(3.8) 



Z{t) 



from which we easily recover eq. ([3J^) at large-n when restricting us to the Gaussian. We cannot 
deduce the macroscopic resolvent G(z, t) from eq. ( |3.8| ) for a generic potential as it still depends on 
the matrix R from the diagonalization of J. To disentangle J and we would need more identities 
of the type eq. (p. 7]). However, eq. ( |3.8] ) will be sufficient to obtain the macroscopic resolvent in the 
weak non-Hermiticity limit. It is defined by taking r — * 1, as in the Hermitian limit, but with keeping 
the following product fixed: 

N(l - r 2 ) = a 2 . (3.9) 

Whereas the appropriately rescaled microscopic quantities (see definition ( |3.16| )) will depend on com- 
plex scaling variables and on the parameter a, the macroscopic quantities are projected onto the 
real axis because of r — ► 1. The smoothed resolvent G(z;t) or the macroscopic spectral density, 
p(z) = lmi/v_>oo^~ Rn(z)i only depend on x = Kez since macroscopically Qmz vanishes. This can 
be explicitly seen from eq. ( |3.5| ) where no a dependence occurs. Furthermore, from the defining 
relation between the microscopic and macroscopic density at the origin we have 



lim 



PS(0 = p(0) 



(3.10) 



Taking the results for a Gaussian potential from p5| which depend on £, a and the masses, it can 
be explicitly seen that the macroscopic density p(0) at weak non-Hermiticity is only non-vanishing 
when taking the limit |£| — > oo along the real line R. Therefore the support of the macroscopic density 
collapses from a two dimensional set to an interval, e.g. a = [—2, 2] in the Gaussian case (see eq. 

oo along 



Eq. ( |3.1C| ) also holds at strong non-Hermiticity. There p(0) is non-vanishing when taking the limit |£ 
all directions, indicating that the support extends into the complex plane. 
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(ft.5|)). For the macroscopic quantities on the real line the relation ( 2.16| ) is no longer valid. From the 
loop equation (3.8) at r = 1 and large-n we haveQ 



G(x±ie,t) = —V'{x)^mp t (x), xEaCR, (3.11) 



where 



2t 

d 



V & s Ef^' 92 = 1, (3.12) 

k=l 

p t (x) = £( 2 , fc V(^ ( 3 - 13 ) 

j=l k=0 ^ ' 

All quantities are now real valued functions depending on x = ?R.ez only, with the above eqs. being 
familiar from the Hermitian one- matrix model (e.g. in ||). The quantity r(t) is the limiting function 
of the recursion coefficients r n at large-n in the three step recursion relation generalizing eq. ( |3.4|) . It 
is related to the endpoint c of the eigenvalue support, a = [— c, c] , through c 2 = 4r(i). The determining 
equation following from the normalization condition f a dxpt(x) = 1 is reading 

i£^ffW)' = l. (3.14) 



2 ^ t \ k 

k=l v ' 

With all ingredients supplied we can finally take the microscopic scaling limit of the asymptotic 
differential equation ([T^) at the origin. It is defined as 

Nz = NQRez + iSmz) = £ . (3.15) 

In this limit the complex eigenvalues z approach zero while £ is kept fixed. The rescaling of the 
microscopic correlators is given by 

Ps(£l,...,&) = lim N-^RNiN- 1 ^,...,^ 1 ^) . (3.16) 

We take the smoothed eq. (3.3) only depending on the macroscopic resolvent G(z,t). The latter is 
projected to the real axis, due to the weak non-Hermiticity limit eq. ( |3.9[ ), At the same time we 
perform the microscopic limit (3.15), rescaling the differential equation by A -2 to obtain 

d 2 P t (Z) = t 2 G(0±ie;t) 2 P t (O . (3.17) 

Here we have already removed the e from the argument of the polynomials since they are analytic. 
The subscript t indicates that they explicitly depend on t = € (0, 1] which is kept finite at large-A. 
Taking the sum of eq. ( p. 17 ) for both signs and using eq. ( |3. 11| ) we finally arrive at 

d 2 t Pm = - t 2 7T 2 p t (0) 2 P t (0 , (3.18) 

where the contribution from the potential ^'(0) vanishes. The following form for the asymptotic 
polynomials thus holds: 

*a p n (z = i) "> p t(0 = I f ^ C °T P t ] S n > * = I > (3-19) 

n,N-^oo; t->1 V N J vs/ 1 g(t) sm(tir p t (0)£) n odd ' A v ' 



7 Since the resolvent is no longer denned on the support a of real eigenvalues we add a small imaginary part He. 
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where we have used their parity to select the appropriate solution. The normalization constants f(t) 
and g(t) still have to be determined. Our universal parameter pt(0) is the i-dependent macroscopic 
spectral density for the partition function eq. (|2.1l| ) with potential V(x)/t from eq. (|3.12| ). At t = 1 
it coincides with the macroscopic density p(0) of our model eq. Q2.1| ) in the Hermitian limit (times a 
delta- function in 9mz). We now introduce the following abbreviation, 

7rtp t (0) = u(t) , (3.20) 

identifying the function u(t) = ds/{2\J r(s)) as it occurs in the universality proof of ref. ||. The 
asymptotics of orthogonal polynomials in eq. ( |3l9l ) is thus the same as that of the polynomials on the 
real line ||, replacing the argument by a complex scaling variable £. However, the unknown coefficients 
f(t) and g(t) may still depend on the potential V{x) in a complicated fashion. In particular, they also 
depend on the weak non-Hermiticity parameter a which is not present for real eigenvalues. Another 
difference is that because of the lack of a Christoffel-Darboux formula the correlation functions from 
eq. (2.10) do not only depend on Pt(£) at t = 1 but on an integral over t [24, |25|]. It is therefore a 



non-trivial task to show that the correlation functions of a Gaussian and our potential eq. (|2.5D agree. 



4 Universality of all correlation functions 



As a first step we calculate the prefactors f(t) and g(t) of the asymptotic polynomials eq. ( 3.19 ) 
by taking the microscopic limit ( 3.15| ) of the orthogonality relation ( |2.8[ ). By keeping t = and 
t' = fixed at large- N the Kronecker S^i becomes a delta-function. On the left hand side we have to 
introduce microscopic variables replacing z = jj. In the potential eq. fl2.5|) terms of higher order than 
quadratic are suppressed in the limits Q3.9| ) and ( |3.15 ), leading to exp[— ^(9m() 2 ]. We arrive at 



dMe£dSm.£ e"^ (3m ° 2 P t {£)Pv{?) 



N 



5(t - 1') . 



(4.1) 



Inserting the asymptotic form of the polynomials ( 3.19 ) we first perform the integral over the real 
part, 

Wot I cos(u(t)£)cos(«(£')£*) 



5(u(t)-u(t')) + £(u(t) + u(0)kcosh(2u(t)9mf) • (4.2) 



The second term can be dropped because of the positivity of the spectral density and thus of u(t). 
The equation for the prefactors then reads 



1 7T 



7 5(t - 



l/(*)| 2 
W)\ 2 



d^sva.^ e 



N 2 u(t) 

where the derivative u(t)' occurs. It leads to the following result 



^ : ' s,M ° 2 cosh(2u(i)9m£) = ^5{t-t'), 



(4.3) 



f{t) = g(t) = y/2Nu{t)'a- X 'K- 



(4.4) 



The constant phase of the two functions which is arbitrary has been set to unity. The final result for 
the asymptotics of the polynomials is reading 



lim 

n,N— >oo; t—*1 



Pn I Z 



1 

N 



V2Nu(t)'a- 1 



3 

7T 2 



5 e _^ u ( t) 2 J cos(u(t)0 n even 



sin(n(t)0 n odd ' * N ' ^ 
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The evaluation of the microscopic kernel defined as Ks(£i, £|) = li m 7V-^oo N 2 Kn(£i/N ", £|/iV) is now 
straight forward: 



Ksfa,® = e-^« 3m ^) 2+ ^) 2 )A /■ 1 ^n(t)'e-^ w W 2 [cosKf)a)cos(u( ^ 2 * ) 



avr Jo V2tt 

+ sm(u(t)Zi)sm(u(t)&)} 

= e ^ u ' — / — = e 2 M cosfuf^i — £>)) . (4.6) 

avr Jo V 2vr 

Here, we have replaced the sum in eq. (U) by an integral, E/Io 1 = ^ Jo 1 dt and 

we have substituted 

t — > u(t). From eq. Q3.20| ) u(0) vanishes at the boundary. Our result eq. ( |4.6| ) exactly coincides 
with the microscopic kernel calculated in for a Gaussian potential and is thus universal. The 
known universal sine-kernel || is reproduced in the limit a — ► 0. The universal parameter is the 
macroscopic density of the Hermitian model p(0) at the origin and it resumes the dependence on all 
the coupling constants in our potential eq. fl2.5|) via eq. ( |3.13| ) at t = 1. The derivation presented 
here is an alternative to the one given in |28| for a Gaussian only, where an integral representation of 
the Hermite polynomials was used. A generalization of such a representation is not easy to obtain for 
our general potential considered. 

The universality of the microscopic correlation functions at weak non-Hermiticity follows from eq. 
( [4 .6j ) by simply inserting eq. (2.10) into the definition ( |3.16| ). In taking the microscopic large- N limit 



of eq. (2/7) it trivially carries over to the massive correlation functions as well, as it has already been 
pointed out in p5[ . For explicit expressions of the massive correlation functions we also refer to pq| . 

In [^J] the microscopic correlation functions at weak non-Hermiticity away from the origin have 
also been calculated for a Gaussian. It is not easy to translate our analysis away from the origin. 
Already in the Gaussian case the leading order coefficient of the asymptotic expansion of the Hermite 



polynomials is not sufficient to reproduce [24]. Since already in the Gaussian case we need the full 
finite-n differential equation for the asymptotics of the Hermite polynomials the question of universality 
in the weak non-Hermiticity limit away from the origin is beyond the scope of this article. 

Let us add two more remarks. As we have mentioned in Section Q the whole formalism of orthogonal 
polynomials also holds in the presence of mass terms in the weight function. A natural question is 
why we could not directly derive the orthogonal polynomials including masses from the very same 



differential equation (|3.3| ). The reason is that the mass terms in eq. (2.1) are given as a product of 



determinants which at large- N would also factorize and thus completely drop out of the representation 



eq. ( 2.12 ). However, we know from the real case [13| that the massive and massless orthogonal 
polynomials are different. Therefore our leading order analysis in eq. ( |3.3| ) does not suffice in that 
case and we have circumvented this problem by the use of relation ([^7]). Our second remark concerns 
the universality of the chiral complex matrix model containing Laguerre polynomials in the complex 



plane, as very recently introduced in 27]. The modification inside the partition function eq. ( ]2.1| ) 
to obtain these orthogonal polynomials can be simply achieved by adding a term —(2a + 1) In \z\/N 
to the potential in eq. ( |2.8D (and by squaring the arguments of the Vandermonde) . However, going 
through the derivation we find that in eq. ( |3.3| ) this additional term is of the same order 0(1) as those 
which we neglected. Unfortunately we know from the analysis of the real eigenvalue case as reviewed 
in H that terms of this order are important and do contribute. Our method is therefore not easily 
applicable to the chiral case. 
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5 Conclusions 



We have shown that universality in random matrix models holds also for correlation functions of 
eigenvalues in the complex plane. The proof we have presented holds in the limit of weakly non- 
Hermitian matrices where the parameter a measuring the non-Hermiticity times the size of the matrices 
N is kept fixed when N goes to infinity. Furthermore, we have restricted ourselves to the microscopic 
limit at the origin where complex eigenvalues at zero are magnified. This limit is particularly important 
in the application to QCD as the spectral density (of Dirac eigenvalues) directly measures the chiral 
condensate, the order parameter of chiral symmetry breaking. 

We have shown universality by analyzing the asymptotics of orthogonal polynomials in the complex 
plane using factorization in the large- iV limit. This lead us to a universal differential equation exactly 
as in the case of real eigenvalues. The universality of the kernel of polynomials and thus of all 
correlation functions then followed. The inclusion of an arbitrary number of Dirac mass terms was 
done by expressing massive through massless correlators. 

Our results present only a first step towards the universality of all complex eigenvalues correlations. 
This is not only due to the fact that the number of symmetry classes is by far larger than that of 
matrix models with real eigenvalues. The exist also two fundamentally different large-TV limit for 
matrices with complex eigenvalues: the limit of weak and strong non-Hermiticity where we have only 
investigated the former. The limit of weak non-Hermiticity is particularly important as it smoothly 
connects the know classes of real and strongly non-Hermitian eigenvalue correlations, by taking the 
limits a — > and a — > oo respectively. These two classes would otherwise be completely disjoint. In 
fact by taking the limit a —* oo of our universal, weakly non-Hermitian correlators we have a good 
argument in favor of universality at strong non-Hermiticity as well. 

Apart from the two classes of non-Hermiticity different correlations may also be found in other 
regions of the spectrum, in the bulk away from the origin and at the edge. Another question is that of 
macroscopic universality, which holds for all connected two- and higher Appoint correlation functions 
of real eigenvalues. Such finding may have important consequences in phenomena such as wetting in 
two dimensions related to the inverse potential problem for analytic curves in the complex plane. 

However, we consider universality of chiral complex matrix models, as very recently introduced by 
the author, as most urgent. Such models are aimed to describe the microscopic fluctuations of Dirac 
operator eigenvalues in four dimensional QCD which are complex due to the presence of a chemical 
potential. To answer that question we need to know subleading corrections in N to the differential 
equation of polynomials. We leave this task for future plans. 

Acknowledgments: I am indebted to P. Di Francesco, E. Kanzieper and G. Vernizzi for enlightening 
discussions and correspondence. This work was supported by the European network on "Discrete 
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